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AERODYNAMIC  COMPRESSION  OP  A  GAS  STREAM 


Ye.  D.  Nesterov 
( Kazan) 

Let  us  examine  the  interaction  of  two  streams  of  compressible 
fluid  directed  at  an  angle  to  one  another  (Pigs.  1  and  10)  .  The  main 
stream  of  gas  is  compressed  by  an  active  stream,  which  is  supplied 
through  an  annular  slit  b,  located  in  the  region  of  minimum  cross 
section  of  the  main  stream. 

Ratios  of  aerodynamic  compression  are  derived  for  two  cases: 
a)  a  convergent  channel  (Pig.  l)  and  b)  a  Laval  nozzle  (Fig.  10)  .  The 
difference  between  these  two  cases  is  that  in  a  divergent  channel  the 
separation  zone  is  in  communication  with  the  atmosphere,  and  in  a  Laval 
nozzle  it  is  a  closed  region.  This  gives  unique  conditions  for  solu¬ 
tion  of  the  problem. 


Basic  Premises 

1.  Let  us  consider  the  interaction  of  streams  at  ratios  of  total 
pressure  corresponding  to  the  velocity  coefficients  of  the  active 
stream  X3  <  1.0,  and  at  relative  flow  rates  of  the  active  stream 
AG3  =  fp  <  °-°6- 
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2.  Let  us  derive  the  fundamental  ratios  in  the  absence  of  ex¬ 
ternal  heat  transfer  within  the  framework  of  one-dimensional  theory. 

3.  As  in  Martin's  work  [3],  let  us  assume  instantaneous  mixing 
of  the  streams. 

4.  The  active  stream  expands  up  to  the  static  pressure  of  the 
main  stream  at  a  jet  inlet  angle  ^  ^  <P  ^  tt  and  up  to  the  pressure  of 

P  TT 

the  partially  retarded  stream  P3  =  — n{  X  t' ~cos<p)  when  0  ^  9  ^  ^  . 

A  Divergent  Channel 
The  diagram  is  shown  in  Pig.  1. 

Additional  Conditions 

A.  Between  the  wall  and  the  surface  of  the  stream  after  mixing 
there  exists  a  stagnant  zone  with  pressure  equal  to  atmospheric. 


Pig.  1. 


B. 


In  cross  section  2-2,  static  pressure  is  p2  =  p^  at  a 


subcritical  flow  regime^of  the  main  stream,  i.e.,  at  the  set  drops: 

K-i  ,  where  k  = 
k 


Pi  _  1 

7^XId)< 


H 


[H±] 


is  the  ratio  of  specific 


heats,  and  p2  =  p2 
when 


F  2  I  Tc^T 

IfttJ 

k 

\ 1  »o,  tt  (*id)  ^  ]* 


at  a  supercritical  flow  regime. 
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Let  us  set  up  a  system  of  equations  which  describe  the  process. 


1.  The  condition  of  cylindricalness 

F , — F2  +  AFh  or  F2  +  1 .  (1) 


2.  Prom  the  condition  of  expansion  of  the  active  stream,  the 
ratio  of  total  pressures  of  the  streams 


* 

«3  = 


-  /7(h) 

mh)-rn\  •  cos  f) 


(2) 


3.  The  continuity  equation:  Gi  +  AP3  =  Gs , 
1 ,  3,  and  2  denote  the  main  stream,  active  stream, 
after  mixing,  respectively. 

Using  gasdynamic  functions  [l],  we  obtain 

1 4-  A  G  ==  ~  •  — 2-  .  .  Fj 

«i  p\  g(h)'  ’ 

where  $•§  =  ^  is  the  temperature  ratio, 

-Yi SrV*- 


where  the  subscripts 
and  the  stream 


(3) 


g.  the  acceleration  due  to  gravity  (9.81  m/sec2  )  ,  and  R  the  gas 
constant  (kg-m/kg-deg)  . 

In  Martin's  work  [3],  the  hypothesis  is  taken  that  the  active 
stream  expands  up  to  pressure  p*.  In  this  case,  when  the  ratios  of 
total  pressures  of  the  streams  w*  ^  1.0,  calculation  cannot  be 
applicable,  since  according  to  this  hypothesis  the  flow  rate  of  the 
active  stream  is  zero,  but  on  the  basis  of  experimental  data  [2,3]  it 
can  be  said  that  interaction  between  the  streams  exists,  i.e.,  the 
active  stream  compresses  the  main  stream. 

4.  The  momentum  equation  for  the  volume  of  fluid  between  cross 
section  1-1  and  2-2:  AfJ. 

A  +  — ’  Vcos («-¥)-  A  +  Ph  a/b: 
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Q 

where  I  =  —  v  +  pP  is  the  total  momentum  [kg],  and  v  is  the  velocity 
of  the  stream  [m/sec].  Using  the  gasdynamics  functions,  we  obtain 

*i 

/(*2)+  *1  •«»•(« --*)— 

_  P*  _  ’  ( *0 

-mwwH+-±..F9-nv, 

pi 

where  acr  =  ~^/  gRT*  is  the  critical  velocity  [m/sec] . 

5.  The  energy  equation 

^L.A6-»r+i=a(i+Ao:,)j-»2'  (5) 

Ax  Al 


where  A  = 


k+1 

2(  k  —  1] 


2k 

k  +  1 


.R. 


6.  The  condition  for  relationship  between  total  pressures. 

A.  Plow  without  losses:  total  pressure  after  mixing  should  be 
determined  from  the  condition  of  constancy  of  entropy  S  =  const, 
which  leads  to  the  following  equation  at  Ri  =  R3: 


1  4-  AG, 
AG, 


•In 


^ITCIL 


p'M 


In  ■ 


K  f  »» 
_•*! — 1  { 

IT  i 


.  _  * 

From  calculation  made  for  k3  =  ki,  £ 3  =  1.0,  Ag3  =  0.05*  tt3  =  0.75  to 
1.4,  we  have  a  difference  of  +  1.5%  between  p2  and  pi.  Therefore,  let 
us  take  p|  =  pi  for  calculation.  * 

B.  Plow  with  losses  in  the  main  stream.  In  this  case  p^  =  a. 

l 

where  a  is  determined  by  experimental  data  as  a  function  of 

a—  * 

AG3,  7T3,  <p. 

7*  The  ratio  of  flow  rates  of  the  streams 


m, 


V|V,(X»\-.A/y 
#,  •  q  (h) 
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(6) 


where  h3  is  the  flow-rate  factor  of  the  active  stream,  selected  from 
experimental  data.  It  is  a  function  of  the  geometry  of  the  slit  (inlet 
conditions)  and  the  Reynolds  number  Re  -  ,  where  v  is  the  kinematic 

coefficient  of  viscosity  [m2/sec]  and  _b  is  the  width  of  the  slit. 

With  smooth  input  to  the  slit,  |x3  is  close  to  1.0. 

As  a  result,  we  have  a  system  of  six  equations  with  eight  unknowns 

af3  ,  $3,  -jg,  f2,  AFjj  ,  n  ( • 

The  gas  constant  and  the  ratio  of  the  specific  heats  of  the 
mixture  are  determined  by  the  formulas 


*1  p  ■  — 2—  AG.-/?, 
K.  +  4G,/?,  L  *,-1  *,-» 

A?  ]  — -  —  ,  «:  - 


1  +  45, 


'  AG, 


A, 


+  ~k  — 


The  two  additional  equations  are  functions  of  the  energy  aggregate  of 

the  main  stream.  Let  us  consider  the  following  schemes. 

.  * 

1.  Apparatus  operating  by  suction.  Let  n  ( =  const  Ti  =  const 
The  ratio  of  the  flow  rates  of  the  main  stream  in  operation  with 
compression  and  without  it,  at  a  subcritlcal  flow  regime  of  the  main 
stream 

* 

p _  Gx _  q  (*i) 


at  a  supercritical  flow  regime 


where  the  subscript  ( 0)  pertains  to  the  parameters  of  the  main  stream 
in  the  absence  of  aerodynamic  control. 

2.  The  main  stream  is  supplied  by  a  centrifugal  compressor. 

Two  versions  are  possible  here. 
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Pig.  2. 

A.  The  number  of  revolutions  is  held  constants  compression  is 
accomplished  by  the  interaction  of  the  two  streams.  On  the  charac¬ 
teristic  of  the  compressor,  the  point  of  joint  operation  is  shifted 
from  position  0  to  position  B.  Assuming  that  the  characteristic  of  the 
compressor  is  horizontal,  we  obtain  for  the  degree  of  compression  the 
condition 

*  P\ 

**  **=  —7-  =  const, 

Pbx 

where  p£x  is  the  total  pressure  of  the  stream  at  the  input  to  the  com¬ 
pressor. 

Under  the  condition  of  adiabatic  efficiency  of  the  compressor 
q*  =  const,  we  have  T*  =  const,  i.e.,  the  problem  reduces  to  the 
previous  case. 

B.  The  number  of  revolutions  is  varied.  Consequently,  on  the 
characteristic  of  the  compressor,  the  point  of  joint  operation  is 
shifted  forward  ( the  number  of  revolutions  increases)  from  position 

0  to  C,  and  backwards  (the  number  of  revolutions  is  decreased)  to  C", 
instead  of  OA'  and  OA" ,  respectively,  with  choking  unchanged. 

Assuming  efficient  operation  of  the  compressor  proportional 
to  the  square  of  the  number  of  revolutions,  and  q*  =  const,  we  have 
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3.  The  main  stream  is  supplied  by  an  axial  compressor.  In  this 
case  let  us  examine  a  variation. 

The  number  of  revolutions  is  held  constant}  compression  is 
accomplished  by  the  interaction  of  the  two  streams.  On  the  character¬ 
istic  of  the  compressor,  the  point  of  joint  operation  of  the  system 
is  shifted  from  position  0  to  B.  Assuming  that  the  characteristic 
of  the  compressor  is  vertical,  for  =  const  we  can  write 


’'At,-1'*.  <?(X,) 


"V 


_»  90) 

*  »<*.) 


i f  vt  --  »<v>fcfrr- 

j/f.  ~ 


k-l 

IL  —  ”*.  ~(1~  ^ 

•  »-t 

^2  *  ft 

"ft.  -  (*  “  ’!*) 


In  reality,  the  efficiency  of  a  compressor  is  usually  somewhat 
lower.  In  the  presence  of  such  compressor  characteristics,  the  calcu¬ 
lation  for  both  cases  is  carried  out  more  strictly. 

By  analogy  with  the  examples  examined,  we  can  obtain  additional 
conditions  for  the  main  stream  for  the  use  of  aerodynamic  control  on 
direct-flow  or  rocket  engines.  As  an  example,  let  us  examine  the  solu¬ 
tion  for  apparatus  operating  by  suction  ( applicable  to  experiments 
[3])  • 

c 

The  ratio  of  specific  heats  — =  kx  =  k3  =  1.4. 

cv 

The  temperature  ratio  5- 2  =  1.0. 

3 
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The  drop  in  the  main  stream  P*/Pjj  =  1-89. 

In  this  case  the  problem  reduces  to  the  solution  of  a  system 
of  two  equations  with  two  unknowns,  Xx  and  X3: 


n(\A)- 


q  (Xs) 

* 


— /(xi)  +  k  -a Ft-q  (X,)  COS  (w  —  <p), 

„M  *•  — >A  /7(X,)  =  /7(X^. 

W,_  /7(Xl.cosT)./7(X3)  o 1,0,  °=7^. 

The  solution  is  carried  out  by  the  method  of  successive  approximations. 
Letting  XJ  <  1.0,  we  find  X”  of  the  second  approximation.  Let  us  check 
the  convergence  etc.  up  to  the  required  accuracy. 


The  efficiency  of  aerodynamic  compression  is  conveniently 
evaluated  by  the  specific  area  of  compression 

8=  l~F3--100°/o. 

The  calculation  was  made  for  relative  magnitudes  of  the  slits  Af3  = 

=  0.053,  0.04,  0.027*  and  0.007.  The  results  are  given  in  Pigs.  4,  5 
and  6. 

In  Pig.  4  it  can  be  seen  that  at  low  AG3^;  0.05  for  Af3  =  2. 7$  5 
at  9  =  70°,  the  calculation  is  in  good  agreement  with  experimental 
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data  [3]. 

In  order  to  determine  the  effect  of  the  isentropic  exponent  and 
the  ratio  of  the  stream  temperatures,  calculations  were  made  for 

A  Fa  =0,027,  <p=  90°,  v-=  1,89. 

k  =1,4,  *  =  1,33,  *=1,25. 

»,2=1,5,  1,0,  0,5. 

The  results  are  shown  in  Pig.  7* 


values  of  9  and  Af3. 


Pig.  5*  The  function 
for  slit  Af3  =  0.7$. 
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Using  the  experimental  dependence  Go  =  f(AGa)  (Fig.  6)  ,  let  us 
calculate  the  losses  in  total  pressure  of  the  main  stream  (Fig.  9) • 

For  aerodynamic  control  of  the  engine  nozzle,  the  inlet  angle  of  the 
active  stream  must  be,  from  the  condition  of  allowable  losses  in  total 
pressure,  90  to  70°. 


Fig.  7*  Variation  in  specific  area 
of  for  divergent  nozzle  and  Laval 
nozzle  at  Ap3  =  2.7$,  <P  =  90°. 
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Let  us  compare  the  calculation  with  experimental  data  [2]  for 
a  slit  Ap3  =  5.0$,  <p  =  90°.  The  depth  of  penetration  of  the  active 
stream  in  Zenukov's  work  [2]  was  estimated  by  the  magnitude  of  the 
aperture  ring  h9  (Fig.  9),  which  exerts  on  the  main  stream  an  effect 
which  is  equivalent  to  aerodynamic  compression,  i.e.,  the  same  flow 
rate  of  air,  with  the  same  total  pressure  and  temperature.  Transition 
from  the  depth  of  the  aperture  ring  h9  to  the  effective  depth  h  ^ 
is  by  the  formula 

Kt — =  i  —  Ki*(i  -  An). 


where  (i,  is  the  flow- rate  factor  for  the  aperture*. 

Calculation  of  the  effective  depth  of  compression  of  the  main 
stream  is  by  the  formula 


"ef 


F, 

1  +  A03  ’ 


at  ki  =  k3  =  1.4;  $3  =  1.0. 


Fig.  8.  Variation  in  total- 
pressure  recovery  factor  as  a 
function  of  AG3  for  ^F3  = 

=  2.7$,  <P  =  70°. 


Pig.  9-  The  dependence  h  _= 

=  f ( 7 r*)  for  AP3=  5.0$,  <p  =  90° 


*Rules  27  to  54  on  the  use  and  checking  of  flow-rate  meters  with 
normal  apertures,  nozzles  and  Venturi  tubes.  Committee  on  Standards, 
Measures  and  Measuring  Instruments,  Under  the  Council  of  Ministers, 
USSR,  1956. 
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The  result  of  the  calculation  is  shown  in  Pig.  9.  The  agreement 


of  the  calculation  with  the  experiment  is  satisfactory. 

Laval  Nozzle 


The  diagram  is  shown  in  Pig.  10. 

Additional  premises:  a)  the  stagnation  zone  is  not  extended 
beyond  the  cylindrical  section,  and  b)  the  pressure  in  the  stagnation 
zone  is  assumed  equal  to  the  static  pressure  in  cross  section  2-2. 

For  a  Laval  nozzle  without  a  cylindrical  section,  it  is  necessary  to 

have,  from  experiments,  the  pressure  in  the  stagnation  zone  for  various 

_  * 

values  of  Ag3,  9,  and  tt3  . 

Let  us  set  up  a  system  of  equations  for  this  case,  assuming 

ki  =  k3. 

1.  The  equation  of  continuity  between  cross  sections 

a)  1  —  min  q(Xj  =  i/lh 

(1  +  ao3)  V  t\ 

P* 

where  cr  =  — —  is  the  total-pressure  recovery  factor; 

*1  *  — 

P  W 

b)  min  —  2  a2  =  tS  =  min  , 

Pmin  qTO 

where  1  —  a2  are  the  losses  in  total  pressure  in  the  section  min  —  2. 

2.  The  momentum  equation  for  the  volume  of  fluid  between  the 
cross  sections 

a)  1  —  min  /i  +  ~  v3  ’  cos  (*  —  f).®8  lm'n  Jrl>3 '•  ' 8  ’ 

Using  gasdynamic  functions,  we  obtain 

z  (X#)  +  -i- .  AO,X,  •  6,  •  COS  («  -  f )  =  ( 1  +  AO,)  #2  •  5  (^2); 
k+  1 

b)  min  -  2  I2  =  Ijnin  +  P3.3  *  AF3>3. 

Taking  point  1  (b)  into  account,  we  obtain 


4?,.a=i-p...-=2[*-Trr  1- 
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active  stream 


-Ur*l\rL. 


n  zone 


Fig.  10 


3.  The  energy  equation,  the  ratio  of  the  flow  rates  of  the 
streams,  and  the  condition  of  expansion  of  the  active  stream  are 
analogous  to  formulas  (5),  (6),  and  (3),  respectively. 

As  a  result  we  have  a  system  of  seven  equations  with  seven 
unknowns:  \lt  X3,  X2,  Af3>3  ,  Fmin,  a2,  and  $2.  The  parameters 
Ag3,  7 r*,  cp,  and  &3  are  given  by  the  conditions  of  the  problem,  and 
a  is  determined  from  experimental  data. 

Therefore,  all  parameters  are  determined  uniquely.  Calculation 
is  by  successive  approximations,  similarly  to  the  problem  examined 
above . 

The  graph  in  Fig.  7  gives  the  results  of  the  calculation  of  the 
specific  area  of  compression  6  as  a  function  of  the  ratio  of  total 
pressures  for  Af3  =  0.027.»  9  =  90°,  ki  =  k3  =  k  =  1.33,  ^3  =  1.0. 
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